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Abstract. We develop a general approach for studying the cumulative probability distribution function of 
localized objects (particles) whose dynamics is governed by the first-order Langevin equation driven by 
superheavy-tailed noise. Solving the corresponding Fokker-Planck equation, we show that due to this noise 
the distribution function can be divided into two different parts describing the surviving and absorbing 
states of particles. These states and the role of superheavy-tailed noise are studied in detail using the 
theory of slowly varying functions. 
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1 Introduction 

The Langevin approach, i.e., incorporation of the noise 
terms into deterministic equations of motion to describe 
the stochastic dynamics of particles, is one of the most ef- 
fective tools for studying the effects induced by fluctuating 
environment. If these terms arise from the time derivative 
of the noise generating processes, i.e., random processes 
with independent and identically distributed increments, 
then the solution x(t) of the corresponding Langevin equa- 
tion is a Markov process [TJ[2] ■ The main feature of Markov 
processes is that their future behavior depends only on the 
present state, or in other words, the correlation time of 
Markov processes equals zero. Although physical processes 
are usually characterized by nonzero correlation times, the 
Markov processes can be considered as a first approxima- 
tion of physical ones if the correlation time is the smallest 
time scale of the system. 

In the Markovian case, the statistical properties of 
the particle position x(t) are controlled by the distribu- 
tion of increments of the noise generating process. If this 
distribution is described by a scaled probability density 
p(y) (see the next section) then, depending on the asymp- 
totic behavior of the tails of p(y), three qualitatively dif- 
ferent cases can be distinguished. The Gaussian density 
represents the first case in which the process x(t) is con- 
tinuous and the probability density P(x,t) of x{t) satis- 
fies the ordinary Fokker-Planck equation OS] . This equa- 
tion describes a wide variety of noise phenomena, includ- 
ing particle diffusion in external potentials [3111], stochas- 
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tic resonance [5], noise- induced transport [SHZIIB], noise- 
enhanced stability [9l[T0] , cross-correlation effects [TT1IT2] , 
and many others. The Levy stable densities represent the 
second case in which the probability density P(x, t) sat- 
isfies the fractional Fokker-Planck equation [LHIhUfPolfTfJl 
H71ITS] that describes a special class of discontinuous pro- 
cesses x{t), the so-called Levy flights [Il?ll20ll21"ll22] . It has 
been recently shown [53] that the ordinary Fokker-Planck 
equation holds for all probability densities p(y) with finite 
variance, and the fractional one for all heavy-tailed p(y), 
i.e., probability densities with power-law tails and infinite 
second moment . 

Finally, the third case corresponds to the superheavy- 
tailed p(y), i.e., probability densities whose all fractional 
moments are infinite. Because of this unusual feature, the 
distributions with superheavy tails are rarely used in phy- 
sics. Nevertheless, the examples of ultraslow diffusion [2H 
25,26,27 demonstrate the utility of these distributions in 
modeling the systems with an extremely anomalous be- 
havior. Very recently, the usefulness of superheavy-tailed 
distributions has also been demonstrated for the Langevin 
approach [35]. Specifically, using the generalized Fokker- 
Planck equation [T5J , it has been shown that superheavy- 
tailed noise, i.e., noise arising from a generating process 
whose independent increments are distributed with super- 
heavy tails, induces two probabilistic states of a particle, 
surviving and absorbing. Due to this feature, the Langevin 
equation driven by superheavy-tailed noise can be used 
to describe some randomly interrupted processes. From a 
physical point of view, the interruption can be associated 
with the transition of a particle to a qualitatively new 
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state. In particular, such a situation occurs when the par- 
ticle dynamics is accompanied by the absorption of these 
particles. It should be stressed, however, that the Langevin 
equation driven by superheavy-tailed noise and the corre- 
sponding generalized Fokker-Planck equation describe the 
simplest situation when the absorption rate does not de- 
pend on the spatial and temporal variables [28]. We note 
in this context that the particle dynamics in an absorbing 
medium characterized by an arbitrary dependence of the 
absorption rate on these variables can be studied within 
the path-integral approach (see, e.g., Refs. [2U1I5D] ). 

While the main purpose of Ref. [35] was to show both 
analytically and numerically that the surviving and ab- 
sorbing states exist, in this paper we present a detailed 
analytical study of the distribution function of particles 
subjected to superheavy-tailed noise. The paper is orga- 
nized as follows. In Section [5] we describe the model and 
introduce the Langevin equation driven by superheavy- 
tailed noise and the corresponding generalized Fokker- 
Planck equation. The main results are obtained in Section 
[3j Here, using the theory of slowly varying functions, we 
confirm our previous results, derive the probability distri- 
bution function of particles, and establish a general con- 
nection between the noise characteristics and absorbtion 
rate. In Section H] for illustrative purposes, we consider 
a special class of superheavy-tailed densities. Finally, our 
findings are summarized in Section [SJ 



2 Definitions and basic equations 

It has been shown in Ref. [18] that the probability density 
P(x, t) of the solution x{t) of the first-order (overdamped) 
Langevin equation 

dx(t) = f(x(i),t)dt + dr)(i) (1) 

satisfies the generalized Fokker-Planck equation 

Q-P&t) = —^f{x,t)P(x,t)+F- l {P k {t)\nS k ). (2) 

Here, we interpret x(t) [x(0) — 0] as a particle position, 
f{x,t) as a deterministic force and drj(t) as a noise term. 
The noise generating process r](t) [r](0) = 0] is assumed to 
have independent and identically distributed (for a given 
dt) increments drj(t) — rj(t + di) — rj(t). It is also assumed 
that the solution of equation ([2J obeys the initial condi- 
tion P(x, 0) = S(x), where 6(x) is the Dirac S function, and 
the normalization condition J^ D oo dxP(x,t) = 1. The last 
term in @ contains the direct and inverse Fourier trans- 
forms defined as J-{u k } = u k = dx e~ lkx u{x) and 

F^iuk] = u(x) = (1/2tt) J* dke lkx Uk, respectively. 
Specifically, Pk(t) = dxe~ lkx P(x,t) is the character- 
istic function of x(t), i.e., P k (t) = (e~ tkx w), where the 
angular brackets denote averaging over the sample paths 
of the noise generating process r)(t), and S k — (e^ tkn ^) 
is the characteristic function of the random variable r](l), 
i.e., noise generating process at t = 1. 



The main advantages of equation @ are that it (i) is 
valid for all noise generating processes rj(t), (ii) accounts 
for the noise action in a unified way, namely through the 
characteristic function S k of 77(1) , and (iii) reproduces all 
presently known Fokker-Planck equations associated with 
the Langevin equation ([T]) jjJT] • The possible forms of this 
equation, which are determined by the possible forms of 
the characteristic function S k , are restricted by the con- 
dition that the random variable 77(1) is represented as an 
infinite sum of independent and identically distributed in- 
crements At)(Jt) = r)(jT + t) — 7?(jr) when the split time 
t tends to zero, i.e., 

[1/t]-1 

^(l)-lim J2 A ^ T ) (3) 

T ^ 3=0 

([1/t] is the integer part of 1/r). For example, if the prob- 
ability density of the increments Arj(jr) is given by 

p(Ar 1 ,r) = ^p(^L), (4) 
a(r) \a(T)J 

where the probability density p(y) satisfies the condition 
lim e _).o p{y I e) / £ = S(y) and the scale function approaches 
zero as t — > 0, then, depending on the asymptotic behav- 
ior of p(y) at \y\ — > 00, the ordinary or generalized central 
limit theorem can be applied to determine S k [T51IHT] (see 
also Sect. I3.2.1[) . It has been shown in particular that if 
p(y) has heavy tails then S k is the characteristic function 
of Levy stable distributions and, as a consequence, equa- 
tion ([2]) reduces to the fractional Fokker-Planck equation. 
On the other hand, since [TB] 

In S k = lirn - \p ka r r \ - 1] (5) 

T-S-0 T 

with p ka ( T ) = dy e~ lka ^ y p{y), the parameters of the 
last equation can be expressed through the asymptotic 
behavior of the probability density p(y) at \y\ — > 00 and 
the scale function a(r) at r — >• 23 . 

In this paper we focus on the study of the cumula- 
tive distribution function F(x, t) — J_ dxP{x, t) of par- 
ticles subjected to superheavy-tailed noises. These noises 
arise from the noise generating processes characterized by 
the superheavy-tailed densities p(y) whose fractional mo- 
ments M x — J_ oo dy\y\ x p(y) are infinite for all x > 
(the normalization of p(y) implies that Mq = 1). Next 
we restrict ourselves to the symmetric superheavy-tailed 
densities with the following asymptotic behavior: 

p(y) - -h(y) (y -> 00), (6) 
V 

where a positive function h(y) is assumed to be slowly 
varying at infinity, i.e., h{ny) ~ h(y) (y — > 00) for all 
fi > 0. Since for all slowly varying functions the condi- 
tion y x h{y) — > 00 holds as y — > 00 [32], the probabil- 
ity densities characterized by the asymptotic behavior © 
are indeed superheavy-tailed, i.e., M x> q = 00. It should 
be noted, however, that the slowly varying functions h(y) 
are not arbitrary because the condition Mq = 1 implies 
that h(y) = o(l/lny) as y — > 00. 
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3 General results 

3.1 Characteristic function Sk 

For finding the characteristic function Sk of the random 
variable 77(1) we use the basic relation (|S|), which with the 
symmetry property p(—y) = p{y) and the normalization 
condition 2 J Q dyp{y) = 1 can be rewritten as 



function Y(X). In general, according to (|5J), this behavior 
is controlled by a given probability density p{y). However, 
because Y(X) varies slowly, there is a possibility to write 
Y(X) at A — > 00 in a quite general form. Such a possibility 
is provided by the Karamata representation theorem |32) , 
which states that for some I > and all A satisfying the 
condition A > I every slowly varying function L(X) can be 
written in the form 



lnS fe = -21imiy(— 

t^ot \\k\a{T) 



(7) 



1(A) = //(A),-q»( J du^- 



(13) 



where 

rOO 

Y(X) = / dy[l - cos(y/X)]p(y) (8) 
Jo 

is a positive function of A that approaches zero as A — > 
00. In those cases when the indefinite integral Jdyp(y) is 
known (see, e.g., the example in Section[4]) it is convenient 
to use an alternative representation of Y(X): 



Y(X) = I dx sinx / dyp(y). 

J0 J\x 



(9) 



Since the function h(y) in the asymptotic formula ([B]) 
is slowly varying, Y(X) also varies slowly. Indeed, using 
(0, one can easily verify that 



Y(fj,X)=fj,X / dx(l — cos x)p(p,Xx) 
Jo 



1 — POS T 

dx —-h(pXx) ~ y(A) (10) 

x 



(A ->■ 00). Therefore, if k ^ then F(l/|fc|a(r)) in (0) can 
be replaced by F(l/a(r)) yielding lnSfc = — q, where 



q = 2\\m-Y(^ 

T^O T \ a ( T ) 



(11) 



is a non-negative parameter that does not depend on k 
(its physical meaning will be discussed in Sect. 13 . 3[) . In 
contrast, if k — then po = J_ oo dyp(y) — 1 and so In So = 
0. Combining these results and introducing the Kronecker 
delta Sko, we obtain [25] 



lnSfc = -q(l - 5, 



k-ol 



(12) 



Thus, for all symmetric probability densities p{y) with 
superheavy tails the characteristic function of 77(1) has the 
form S k = e - q{1 - Sk0 \ It should be noted that although S k 
differs from e~ q only in one point k — 0, the Kronecker 
delta cannot be neglected because it provides a correct 
normalization of the distribution of 77(1) (see Sect. I3.2|) . 



3.1.1 Representations of q and a(r) 

As it follows from (fTTj) . the parameter q, which accounts 
for the influence of superheavy-tailed noise on the system, 
depends on the asymptotic behavior of the slowly varying 



where g(X) ->j£ (0, 00) and e(A) — » as A — > 00. Since 
the defining property of L(X), i.e., L(pX) ~ L(X), involves 
only the asymptotic behavior of this function, the repre- 
sentation (|13p is essentially non-unique. Using this fact, 
it can be shown [32] that L(X) at A — > 00 can always be 
represented in a simpler form 



L(X) - gexp 



dve(e v ) 



(14) 



'In 2 



Thus, assuming that Jdve(e v ) 
varying function 1^(A) we obtain 



-^(u), for the slowly 



Y(X) ~ cexp[-<?(lnA)] 



(15) 



(A — > 00) with c = g cxp[<P(\nl)]. An explicit form of the 
function <P(v) depends on p(y), but since Y(X) tends to 
zero when A increases, the condition <P(v) — > 00 as v — > 
00 must hold in all cases. Using (fTTj) and (fTS"|) . for the 
parameter q we find the desired general representation 



2c 

q = — , 

r 



where r is a time scale parameter defined 



as 



r = lim r exp 



<P\ In- 



a(r) 



(16) 



(17) 



Depending on how fast the function a(r) approaches zero 
as t — y 0, the parameter r can take the values from the 
whole interval [0, 00]. If a(-r) tends to zero so rapidly that 
r = 00 then such noise does not affect the system at all 
(q = 0). In the opposite case, when r = 0, the noise action 
is so strong that the system reaches the final state at t = 
+ , i.e., immediately (q = 00). 

In what follows we are interested in a non-trivial action 
of superheavy-tailed noise. This case is specified by the 
condition < r < 00 and occurs only if the scale function 
a(r) has a proper asymptotic behavior at t — > 0. In order 
to find this behavior, we assume that for large enough 
ip there exists the inverse function v — <P^ 1 (ip) of the 
function cp — <I>(v) (since <P(v) — > 00 as v — > 00, in this 
case •P^ 1 (ip) -> 00 as -> 00). Then the limit in (fl7|) is 
equal to a given value r if the asymptotic behavior of a(r) 
at t — y is given by 



a(r) ~ exp 



-^-Mln 



(18) 
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We note that the asymptotic formula (fTS)) is not unique 
because, due to the presence of the logarithmic function 
in ([TT)) , any scale function of the form a(r)6(r) with b(r) 
satisfying the condition 



t exp \ <P 



{<Z> (p-^ln-) -ln6(r) } 



(19) 



(r — > 0) also leads to the same parameter r. It is important 
to emphasize, however, that this lack of uniqueness is of 
no importance because the parameters c and r completely 
characterize the noise action. Therefore, at fixed c and 
r the Langevin equation ([TJ driven by superheavy-tailed 
noise and the corresponding Fokker-Planck equation (T2J| 
are well defined and the function b(r) can always be chosen 
as b(r) = 1 without loss of generality. 



3.2 Distribution function of »7(1) 

In general, using (T2J and (|12[). it is possible to establish 
the main features of the probability density P(x, t) aris- 
ing from the action of superheavy-tailed noise [28] . How- 
ever, because the random variable 77(1) is responsible for 
all these features, it is reasonable to study the statistical 
properties of 77(1) in more detail. This will also permit 
us to obtain the probability distribution function of the 
particle position, F(x,t), in a consistent way. 

In order to gain more insight into the properties of the 
random variable 77(1), we calculate its distribution func- 
tion F(rj). To this end it is convenient to temporarily re- 
place the Kronecker delta 5ko in the characteristic function 
Sk = e~ q { x ~ 5h0 ' of ?y(l) by a symmetric and smooth func- 
tion of k, Akin), which depends on a positive parameter 
re in such a way that lim K ^o Ak{n) = Sko- We assume that 
Sk(n) = e~ q ^~ Ak ^ is the characteristic function of the 
random variable 77(1; re) whose probability density S(r); k) 
and the corresponding distribution function F{ji; re) are 
given by S(tj; re) = (1/2tt) JT^ dk e ik ^S k {n) and ^(77; re) = 

dr]S{r}; re), respectively. If the distribution function 
^(77; re) is known then the desired distribution function 
can be determined as F(rj) = lim^^o ^(77; re). 

The choice of the function /Afc(re) is, of course, not 
unique. For example, it can be chosen as 



^fc(re) = - In 

q 

(sech x = 1/ cosh x) or 
1 



l + (e 9 -l)sech( — 
1 2re 



A k M = - In 1 + (e« 

q 



l)e 



-\k\/ K 



(20) 



(21) 



However, in the limit re — > the explicit form of Zifc(re) 
is of no importance. Therefore, without loss of generality, 
we consider Z\/c(re) from (|2"0")) for which the characteristic 
function of 77(1; re) has the form 



S k {K) 



(1 



> q ) sechl — 
1 2re 



(22) 



The probability density and the distribution function that 
correspond to this characteristic function are given by 



5(77; re) = e- q S(T]) + (1 - e~ q ) - sech rery (23) 

7T 



and 



F(r r , re) = e- q F d (f]) + (1 - e- q )F( K r)), (24) 
respectively. Here, 



F d (v) 



0, 77 < 

1, 77 >0 



(25) 



is the step function describing the degenerate distribution 
localized in the point 77 = 0, and 



F(kti) = - arctan(e K?? ) 



(26) 



is the distribution function which describes the hyperbolic 
secant distribution. 

According to (|24l) , the distribution of the random vari- 
able 77(1; re) is the mixing of the degenerate and hyper- 
bolic secant distributions. In contrast to the former, which 
does not depend on the parameter re, the latter strongly 
depends on this parameter. In particular, the maximum 
height of the distribution, dF{Kri)/drf\ n =Q, and the root 
square of the variance, which characterizes the width of 
this distribution, are equal to re/7r and 7r/(2re), respec- 
tively It is clear, therefore, that the probability density 
dF(K,rf)/dr) of the hyperbolic secant distribution tends to 
zero as re — > 0. On the other hand, as it follows from 
(f2"u) . this occurs in such a way that F(— 00) = and 
F(oo) = 1 for all re > 0. Assuming that at re — > these 
conditions hold as well, for the limiting distribution func- 
tion F\(rj) = limK^o F(kt]) we obtain 



0, 77 = —00 
F,(, 7 ) = { 1/2, M <co 

1, 77 = 00. 



(27) 



This particular distribution function describes a normal- 
ized random variable whose probability density equals zero 
for all |?/| < 00. From a formal point of view, it can also 
be considered as a discrete variable that takes only two 
values, 77 = —00 and 77 = 00, with probability 1/2 each. 
Thus, the probability distribution of 77(1) characterized 
by the distribution function ^(77) = lim K _i.o F(r]; re) is the 
mixing of the degenerate distribution (|25[) taken with the 
weight e~ q and the limiting distribution (|2"T|) taken with 
the weight 1 — e~ 9 , i.e., 



^(77) = e- q F d (r)) + (1 - e- q )F x {rj). 



3.2.1 Connection with limit theorems 



(28) 



Representing the increments Arjljr) of the noise generat- 
ing process 77 (t) as Ar)(jT) = yj/c n with c n = l/a(r) and 
n = [1/r] + 1, from ([3]) one obtains 



hm y 

n— >oo — 4 C n 



(29) 
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where yj are the independent random variables distributed 
with the same probability density p(y). We note that in 
probability theory the partial sums of a more general form, 
i.e., Y^j=i Vj I ' c n~ d n , is usually considered. But in our sit- 
uation d„ = because a class of probability densities p{y) 
is restricted by the condition \\to. € ^q p{y / e) / e = 5(y). If 
the variance of p{y) is finite then, according to the cen- 
tral limit theorem |33j . c„ oc n 1 / 2 and the distribution of 
77(1) is Gaussian. In this case air) oc r 1 / 2 and for all such 
p(y) the generalized Fokker-Planck equation @ reduces 
to the ordinary one |23j . In contrast, if p{y) with heavy 
tails belongs to the domain of normal attraction of a given 
stable distribution characterized by an index of stability 
a £ (0, 2) then, as the generalized central limit theorem 
suggests [33 , n oc n 1 /" and 77(1) is described by this sta- 
ble distribution. It has been shown [33] that in this case 
a(r) oc T X / a [in special situations, there also exist two dif- 
ferent forms of the scale function a(r)] and equation ([2]) 
reduces to the fractional Fokker-Planck equation. 

Finally, the above derived results show that for all sym- 
metric superheavy-tailed densities p(y) characterized by 
the asymptotic behavior ([5]) the distribution function of 
the random variable 77(1) has the form (|28[). It is important 
to emphasize, however, that in contrast to densities with 
finite variance and densities with heavy tails there is no 
universal scale function a(r) for superheavy-tailed p{y). 
On the contrary, as is clear from ([§]), (fT5"|) and (fT5)) , in this 
case different p{y) lead in general to different a(r). This 
is in accordance with the well-known fact |33) that only 
stable distributions have domains of attraction. Thus, the 
result (j2"5)l can be interpreted as follows: For every sym- 
metric superheavy-tailed p{y) the scale function a(r) can 
always be chosen so that the distribution function of 77(f) 
is given by (|2"5)l with q £ (0, 00). 



3.3 Solution of the generalized Fokker-Planck equation 

Our next step is to solve the generalized Fokker-Planck 
equation ^ for the characteristic function S k = e~ q ( x ~ 5h ° > 
that corresponds to superheavy-tailed noise. To this end 
we first solve this equation for the characteristic function 
(|2"2"|) and then take the limit n — > 0. In order to clarify our 
approach, let us designate the solution of equation ([2J in 
which S k is replaced by Skin) as P(x,t;n). Then we de- 
compose it as P(x,t; k) = V(x,t) + A(x,t; k), where the 
terms V(x, t) and A(x, t; n) are governed by the equations 



■^f(x,t)V(x,t)-qV(x,t), 



(30) 



dt 



A(x, t; k) 



Ox 



f(x, t)A(x, t; k) — qA(x, t; k) 



+ qT- 1 {[P k {t)+A k (t;K)]A k (K)} (31) 

and satisfy the initial conditions V(x, 0) = S(x) and 
A(x, 0; k) = 0. We also assume that the normalization con- 
dition for P(x,t; k), i.e., P k (t; n)\ k =o = Vo(t) + Ao(t; k) = 
1, holds for all t and k. Then, with these assumptions, the 



solution of equation (0) can be represented as follows: 

P(x,t)=V(x,t) + MmA(x,t;K). (32) 

K— >0 

An explicit form of the first term in the right-hand 
side of (|32|) can easily be determined from equation ((301) . 
Indeed, assuming that V(x, t) = e~ qt lZ(x, t), this equation 
reduces to 



d_ 

01: 



K(x,t) 



d_ 

dx 



f(x,t)K(x,t). 



(33) 



Its general solution is written as lZ(x,t) = <P[x — z(t)] 1 
where <P(x) is an arbitrary function and z(t) is the solution 
of the equation z(t) = f(z(t),t) (we assume that z(0) = 
and \z(t)\ < 00 for all finite t). Using the initial condition 
1Z(x, 0) = S(x), one finds <P(x) = S(x) and so 



V{x,t) = e~ qt 6[x - z(t)]. 



(34) 



As it follows from this result, the normalization condition 
for V(x,t), V (t) = J^dxVix^), reads V (t) = e~ qt . 

According to ((3T|) and (p4|) . the term A(x 7 t; k) obeys 
a closed integro-differential equation 



d 

A(x, t; k) — ——f{x, t)A(x, t; k) — qA(x, t; k) 



0_ 

Of 



+ q dyA(x - y,n)A(y,t;n) 

J — OO 

+ qe- qt A(x- z(t)-n), (35) 

where A(x;k) = (I/2tt) dk e ikx A k {n) is the kernel 
function. Using the fact that f_ dxA(x, n) — 1, which 
follows from the condition Aq(k) — 1 [see l|20p]. and the 
integral representation of the 8 function S(x) = (1/2tt) x 
J_ dke lkx , from equation (|35[) we obtain a simple equa- 
tion dAo(t; K)/dt — qe~ qt for Ao(t; n) = dxA(x, t; k). 
The solution of this equation satisfying the initial condi- 
tion Aq(0;k) = yields A {t;K) = 1 — e~ qt . Therefore, 
like P(x, t), the normalization condition of A(x, t\ k), i.e., 
Ao(t; k) = 1 — e~ qt , does not depend on the parameter k. 

Because of its complexity, we are not able to solve 
equation (|35l) for an arbitrary k. However, in the case of 
our interest, when k approaches zero, the dependence of 
the solution of this equation on k can easily be found. The 
key point is that the kernel function A(x, n) at \x\n <C 1 
can be approximated as A(x; k) = £,k with 



2 f°° 

£=^- dyln[l + (e 9 -l)sechyl. 
77 Q Jo 



(36) 



In this approximation, the integral term in (f3"5"|) equals 
q£(l — e~ qt )n and A(x, t; n) at \x — z(t)\n <C 1 is governed 
by a linear differential equation 

d d 

— A{x, t; k) = - Q^f( x i t)A{x, t; k) - qA(x, t; k) + q£n. 

(37) 

Its general solution is written as A(x, t; k) — A p (x, t; k) + 
Ag(x, t; k), where A p (x, t; k) oc k is a particular solution of 
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this equation and A g (x, t; k) is the general solution of the 
corresponding homogeneous equation (when k — 0). Since 
A(x, 0; k) = 0, the latter is also proportional to k and so 
A(x, t; k) oc k. Put differently, the term A(x, t; k) tends 
to zero linearly as k — > 0. In particular, if the external 
force is a constant, i.e., f(x,t) — f, then A(x,t;n) — 
£(1 — e~ qt )n at \x — ft\n <C 1. We note also that while 
A(x ) t; k) approaches zero in the limit k — > 0, the spatial 
region of localization of A(x, t; k), which is determined by 
the condition \x — z(t)\ ~ 1/re, tends to infinity keeping 
the normalization condition of A{x, t; n) fixed. 

Thus, according to (l3"2"j) , and the above properties 
of A(x, t; k) at k — y 0, the distribution function F(x, t) = 
j-ao dxP(x, t) (t < oo ) of particles driven by superheavy- 
tailed noise has the form 

F(x,t) = e~ qt F d [x - z(t)} + (1 - e- qt )F x {x). (38) 

It shows that this noise generates two time-dependent 
probabilistic states of each particle. The first one is de- 
scribed by the degenerate distribution function Fd[x— z(t)] 
and is realized with the probability Po(t) = e~ qt . Since 
the particle trajectory z(t) is not influenced by the noise 
[see (l34l) ]. following the terminology of [28] we call this 
state the surviving state. The second state is associated 
with the limiting distribution function Fi(x) and is real- 
ized with the probability Ao(t) = 1 — er qt . In this state 
the probability to find the particle in any finite interval 
equals zero. The transition to this state implies that the 
particle jumps to plus or minus infinity and, in fact, it 
is excluded from future consideration. Therefore, like in 
[28] . we refer to this state as the absorbing state. Because 
of these features, superheavy-tailed noise acts on parti- 
cles as an absorbing medium characterized by the rate 
of absorption q. Hence, the Langevin equation (JTJ) driven 
by this noise and the corresponding generalized Fokker- 
Planck equation @ describe the overdamped motion of a 
particle in an absorbing medium. Moreover, one expects 
that these equations can also be used to describe some 
other processes whose duration is random. 

3.4 Compound noise 

According to the above results, superheavy-tailed noise is 
characterized by two main features which at first glance 
seem to be contradictory. On the one hand, this noise is so 
strong that the transition into the absorbing state occurs 
so that a particle is immediately transferred to infinity. 
But, on the other hand, before this transition the noise 
does not affect the particle position. In other words, before 
absorption which is random in time the particle motion is 
deterministic and is governed by the equation of motion 
z(t) = f(z(t),t) [see (J3U)]- From a mathematical point of 
view, these features result from both superheavy tails of 
the probability density p(y) and appropriate choice of the 
scale function a(r). 

The case when before absorption the particle motion 
is random can also be incorporated into the proposed ap- 
proach. For this purpose we consider the compound noise 



generating process r](t) = f]x(t) +r)2(t), where rji(t) gener- 
ates superheavy-tailed noise, which models an absorbing 
medium, and 772 (t) generates any other noise (e.g., Levy 
stable noise), which causes the random motion of a par- 
ticle. Assuming that these random processes are indepen- 
dent, we obtain Sk — SifcS^fc, where S\k and S 2 k are the 
characteristic functions of 7/1 (1) and 772(1), respectively. 
Then, proceeding as before, we replace Sik by the char- 
acteristic function (|22[) and write the solution of equation 
([2]) at k 7^ as P(x, = V(x,t) + A(x, t; k), where 
V(x, t) is the solution of the equation 

^P(x, t) = -JLf( x , t)T(x, t) - qV(x, t) 

+ T- 1 {V k (t)lnS 2k } (39) 

satisfying the initial condition V(x,0) — 5(x). This so- 
lution can be written in the form P(x,t) = e~ qt W(x,t), 
where the probability density W(x,t) is governed by the 
generalized Fokker-Planck equation 

-WOM) = —^f(x,t)W(x,t) +F- 1 {W k (t)\nS 2k } 

(40) 

and obeys the initial condition W(x, 0) = S(x) and the 
normalization condition Wo(t) = 1. It is also not difficult 
to verify that, like in the previous case, Ao(t;K) = 1 — 
e~ qt and A(x,t; k) — > as n — > 0. Therefore, for all finite 
times the distribution function of particles subjected to 
this compound noise takes the form 

F(x,t)=e- qt f dxW{x,t) + (1 - e^)Fi(^). (41) 
J — 00 

Using the fact that superheavy-tailed noise plays the 
role of an absorbing medium, the above result can be in- 
terpreted as the distribution function of particles whose 
random motion before absorption is described by the fol- 
lowing Langevin equation: dx{t) = f(x{t) 1 t)dt-\-drj2{t). 
The second term in the right-hand side of (|4ip relates to 
the absorbing state and does not depend on the character 
of the particle motion. In contrast, the particle dynam- 
ics strongly influences the first term which describes the 
surviving state. As is clear from (l4Cfl) and f|41 [) . this influ- 
ence is determined by the deterministic force f(x,t) and 
the characteristic function S2k of 772(1). If i] 2 (t) = 0, i.e., 
the particle dynamics before absorption is deterministic, 
then S 2 k — 1, W(x,t) = 8[x — z(t)\ and so the distribu- 
tion function (|4"T1) reduces to (|38D . The random motion 
of particles occurs if the noise generating process rj 2 (t) 
depends on time. As it was discussed in Sect. 13.2.11 in 
this case the possible distributions of 772(1) follow from 
the fact that 772(1) is represented as an infinite sum of 
independent and identically distributed increments. For 
example, if the distribution of the increments Ai] 2 (jT) 
is symmetric and has heavy tails then the random vari- 
able 772(1) is distributed with the symmetric Levy sta- 
ble distribution whose characteristic function is given by 
S 2 k = exp(— 7|fc|"), where a and 7 are the index of stabil- 
ity and scale parameter, respectively. Defining the Riesz 
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derivative as d a h(x)/d\x\ a = -F^^k^ht} [33], in this 
particular case equation (|40)) reduces to the fractional 
Fokker-Planck equation piHil[T51[TCl[r71fT5] 

^W(x, t) = -^f(x, t)W(x, t) + 7^W(z, t), (42) 

which together with the distribution function (l4"Tj) de- 
scribes Levy flights (if < a < 2) and Brownian motion 
of particles (if a = 2) in an absorbing medium. 

Equation (|4"2l can be solved analytically in some sim- 
ple cases. In particular, if /(x, t) is a linear restoring force, 
i.e., f(x,t) — ~bx (b > 0), then the solution of this equa- 
tion is the Levy stable density [T5] 



W{x,t) = J'- 1 {exp(-c 2 (i)|fc| Q )} 



(43) 



characterized by the same index of stability and the time- 
dependent scale parameter c 2 (t) = 7(1 — e~ abt )/ab. At 
a = 2 it reduces to the Gaussian density 



W{x,t) 



20Fc(t) 6XP ( 4c 2 it))' 



(44) 



and the distribution function (|41[) for Brownian particles 
in an absorbing medium takes the form 



F(x,t) 



-it 



■ erfc 



2c(t) 



+ (1 - e-^F^a:), (45) 



where erfc(a;) = (2/t/tt) J 00 dye y is the complementary 
error function. 



4 Particular class of superheavy-tailed 
distributions 

As is shown in the previous section, the distribution of 
particles subjected to superheavy-tailed noise is the same 
as their distribution in an absorbing medium character- 
ized by the rate of absorption q. According to (fTTj) . this 
rate is determined by both the asymptotic behavior of the 
probability density p{y) at y —> 00 and the asymptotic 
behavior of the scale function a(r) at r — >■ 0. Since p(y) is 
assumed to be given, the noise action is non-trivial, i.e., 
< q < 00, if the asymptotic behavior of a(r) is given 
by (|18[) . It should be noted that finding the scale function 
a(r) at small r is important not only for determining the 
absorption rate q, but also for the numerical simulations 
of the Langevin equation [25J . 

In order to illustrate how our general results can be 
used for the determination of q and a(r) at r — > 0, we 
consider a two-parametric class of superheavy-tailed dis- 
tributions that is described by the probability density 



p{y) 



(y- l)(lns) 1 



2( S + |y|)[ln( S +|y| 



(46) 



Since p(y) is assumed to be normalized and non-negative, 
the parameters v and s must satisfy the conditions v > 1 



and s > 1. According to (|4"6"|) . the probability density p(y) 
is symmetric, has a single maximum located at y — 0, and 
its asymptotic behavior at y — > 00 is given by (O with 



(v- l)(lns) 1 
2(lnj/)" 



(47) 



For the probability density (|46p the inner integral in 
([9]) can easily be calculated yielding 



(Ins)"- 1 f 00 
Y(A) = [ -^— I dx 



smi 



o ' [ln(s + Xx)]"- 1 ' 



(48) 



To find the leading term of the asymptotic expansion of 
Y(X) as A — > 00, we use the identity 



1 - e~ ap f p 

= / due-™ 

a 



(49) 



with a = [ln(s + Ax)]" -1 and p — (lnA) 2_I/ . Introducing 
the variable of integration z = uj p and taking into account 
that ap ~ ln(xA) ~ In A if A — > 00 and < x < 00, from 
(gU) we obtain 



1 



1 



[ln(s + Ax)]"- 1 (In A)"- 2 Jo (Ax) 



dz- 



1 



(50) 



Then, substituting (|50")l into (14*51) and using the integral 
formula J °° dxx~ z sinx = r(l — z) cos(7rz/2), where r(z) 
is the gamma function, we arrive to 



Y(X) 



(Ins) 1 



2(lnXy- Jo 



dze 



— z In A 



r(l-*)cos— . (51) 



Since A — > 00, the main contribution to the integral in 
(|5ip comes from a small vicinity of the lower limit of inte- 
gration. Making use of this fact, one can easily verify that 
this integral is asymptotically equal to 1/lnA, and so the 
asymptotic formula (|51[) takes the form 



Y(A) 



1 / In t 



2 V In A 



(52) 



We are now in a position to find the asymptotic rep- 
resentation of the functions <P(v) and ^ -1 (i/?) in the ref- 
erence case [we remind that according to (fT6")) - ([T8l these 
functions determine the absorption rate q and the scale 
function a(r)]. By comparing (l52"j) with (fT5)) . one obtains 
c = (\ns)' / - 1 /2 and $(v) ~ {v- l)ln« (u -> 00). The first 
result leads to the absorption rate 



(In s 



(53) 



and the second permits us to find an explicit asymptotic 
formula for the inverse function: <? -1 (<£>) ~ exp\ip/{v — 1)] 
(tp — > 00). Using this formula, from (fT5)l we get 



a(r) ~ exp 



exp 



-J- in I 
1/ — 1 r 



(54) 
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as t — ¥ 0. Thus, the noise action is specified not only by 
the parameters v and s, which characterize the superheavy- 
tailed probability density (pttjj) . but also by the parameter 
r, which characterizes the asymptotic behavior (jM]) of the 
scale function a(r). We note also that at v = 2 the formu- 
las (I5U)) and ([53]l are reduced to those derived in Ref . [35J , 
i.e., q = (In s)/r and a(r) ~ exp(— r/r). 

We complete our analysis by comparing the sample 
paths of the solutions x(t) of the Langevin equation ([T} 
driven by different noises. If the transition probability den- 
sity of the noise generating process is given by Q then, de- 
pending on the probability density p(y) , the sample paths 
can show three qualitatively different behaviors. Specifi- 
cally, if the variance of p{y) is finite then the sample paths 
of x(t) are random and continuous (e.g., as in the case of 
the Wiener process). In contrast, the sample paths which 
correspond to the probability density p(y) with heavy tails 
are random and discontinuous (e.g., as in the case of the 
Levy process). Finally, if p(y) has superheavy tails then 
the corresponding noise does not influence x(t) up to a 
random time t tr whose average value, t tr = f Q tdAo(t), is 
expressed through the absorption rate as t tr = X/q (here 
dAo(t)/dt = qer qt is the probability density of the sur- 
viving time). At t = t tv an infinite jump of x(t) occurs 
in the positive or negative direction, which is interpreted 
as the transition of a particle to the absorbing state (see 
also Ref. [IS]). Thus, in this case the sample paths are the 
randomly interrupted realizations of the same determinis- 
tic function z(t) [z(0) = 0] satisfying the motion equation 
z(t) = f(z(t),t) (in the case of compound noise the sample 
paths are the randomly interrupted random processes, see 
Sect. I3.4p . It is this property of the sample paths which 
makes it possible to model a number of physical processes 
interrupted at random times by the use of the Langevin 
equation driven by superheavy-tailed noise. 

Examples of these processes can be found, e.g., in sys- 
tems where the particle dynamics displays a qualitative 
change at some random time. One of them, the parti- 
cle propagation in an absorbing medium characterized by 
a constant absorption rate, is studied in this paper. We 
note that because all superheavy-tailed noises generate 
the absorbing state, for the numerical simulations of the 
Langevin equation ([1]) driven by this noise the probability 
density p(y) can always be chosen as in pfjj) . The other ex- 
ample is the position of a fluid particle in a boiling liquid. 
In this case, the interruption of the process corresponds to 
the transition of a particle to the vapor state. One more 
example is the position of a free electron in a semiconduc- 
tor. Here the interruption of the process occurs due to the 
recombination of the electron-hole pair. 



5 Conclusions 

We have developed a general approach for studying the 
statistical properties of particles driven by superheavy- 
tailed noises. These noises arise from the noise generating 
processes whose independent increments are characterized 
by the absence of finite fractional moments. Due to this 



feature, the distribution of particles has been studied in 
detail. Our approach is based on the generalized Fokkcr- 
Planck equation @ which corresponds to the Langevin 
equation ([TJ. Within this framework we have established 
that the distribution of particles is the same as in the case 
of particles moving in an absorbing medium characterized 
by the rate of absorption q. In other words, superheavy- 
tailed noise plays the role of an absorbing medium. 

By analyzing the solution of the generalized Fokker- 
Planck equation, we have shown that the distribution func- 
tion can be divided into two parts which describe two 
probabilistic states of each particle. These states, surviv- 
ing and absorbing, are induced by superheavy-tailed noise 
and are realized with the probabilities e~ qt and 1 — e~ qt , 
respectively. In the surviving state the particle dynamics 
occurs in such a way as if this noise is absent. Therefore, 
the distribution of particles in this state is non-degenerate 
or degenerate depending on that is the driving noise com- 
pound or not. We have derived the corresponding distri- 
bution functions for both these cases. An infinite jump of 
a particle, which occurs at a random time under the action 
of superheavy-tailed noise, is interpreted as the transition 
of this particle to the absorbing state. We have shown 
explicitly that the distribution of particles in this state is 
characterized by the condition that the probability to find 
these particles in any finite interval equals zero. 

Using the theory of slowly varying functions, we have 
derived a general representation for the main character- 
istic of superheavy-tailed noise, i.e., the absorption rate 
q. This representation depends on the asymptotic behav- 
ior of the superheavy-tailed probability density p(y) at 
y — > oo and the scale function a(r) at r — > 0. The asymp- 
totic behavior of a(r) is controlled by the asymptotic be- 
havior of p(y), which is assumed to be known. Finally, as 
an illustration of our results, we have found the absorption 
rate and the scale function for a wide class of superheavy- 
tailed densities. 
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